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ABSTRACT 

We analyze kinematic data of 41 nearby (z < 0.1) relaxed galaxy clusters in terms of 
the projected phase-space density using a phenomenological, fully anisotropic model of 
the distribution function. We apply the Markov Chain Monte Carlo approach to place 
constraints on total mass distribution approximated by the universal NFW profile and 
the profile of the anisotropy of galaxy orbits. We find the normalization of the mean 
mass-concentration relation is c = 6.9^q'y at the virial mass M v = 5 x 1O 14 M . By 
comparison with the calibration from cosmological iV-body simulations it is demon- 
strated that this result is fully consistent with as from WMAP1 data release and 
agrees at ~ ltr level with that from WMAP5. Assuming a one-to-one correspondence 
between eg and the normalization of the mass-concentration relation in the frame- 
work of the concordance model we estimate the normalization of the linear power 
spectrum to be cr 8 = O.Qlj^os- Our constraints on the parameters of the mass profile 
are compared with estimates from X-ray observations and other methods based on 
galaxy kinematics. We also study correlations between the virial mass and different 
mass proxies including the velocity dispersion, the X-ray temperature and the X-ray 
luminosity. We demonstrate that the mass scaling relations with the velocity disper- 
sion and the X-ray temperature are fully consistent with the predictions of the virial 
theorem. 

We show that galaxy orbits are isotropic at the cluster centres (with the mean 
ratio of the radial-to-tangential velocity dispersions a r /ag — 0.97 ± 0.04) and radially 
anisotropic at the virial sphere (with the mean ratio a r /<rg — 1.75+°'^). Although the 
value of the central anisotropy appears to be universal, the anisotropy at the virial 
radius differs between clusters within the range 1 < (cy/cg) < 2. 

Utilizing the Bautz-Morgan morphological classification and information on the 
prominence of a cool core we select two subsamples of galaxy clusters corresponding 
to less and more advanced evolutionary states. It is demonstrated that less evolved 
clusters have shallower mass profiles and their galaxy orbits are more radially biased 
at the virial sphere. This property is consistent with the expected evolution of the 
mass profiles as well as with the observed orbital segregation of late and early type 
galaxies. 

Key words: galaxies: clusters: general - galaxies: kinematics and dynamics - cos- 
mology: dark matter 



1 INTRODUCTION 



Kinematic data for galaxy clusters offer a unique possi- 
bility to study their mass distribution and the orbits of 
the galaxies within them. Depending on the quality of 
data sample and the purpose of the analysis a number of 
methods of data modelling has been propos ed, from the 
mass estimate based on the virial theorem (jHeisler et al.l 



1 19851 : iGirardT et al.|[l998l ) , the Jeans analy s is wit h the two 
even velocity mo ments (lLokas fc Mamonl l2003t ) and the 
caustic technique dDiaferio fc Gellei 1997 ) to the recently 
presented kinematic d eprojection with the mass inversion 
IjMamon fc Boug|2010l ) and the projected phase-space anal- 
ysis with a fully anisotro pic model of the distribution func- 
tion (|Woitak et al.|[2009l ). There is no doubt that the compo- 
nent which is best constrained with all methods of data mod- 
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elling is the mass profile. The analysis of galaxy kinematics 
in nearby galaxy clusters revealed that the mass distribution 
is consistent wi th the universal NFW (INavarro et alJ 1997) 
profile (see e.g. Biviano fc Girardil [20031 ; lLokas etafl l2006t 
IRines fc Diaferiq 20061 ). which was con firmed in other stud- 
ies ba sed on X-ray observati ons (see e.g.lPointecouteau et al.l 
120051 ) or lensing data (e.g. IComerford et al.ll2006h . On the 
other hand, constraints on the parameters of the universal 
density profile of individual clusters from galaxy kinemat- 
ics are available f or a very limited num ber of objects (see 
lLokas et al1l200d : IWoitak fc Lokasll2007l ). One of the main 
aims of this work is to provide robust estimates of these 
parameters for a statistically significant sample of nearby 
(z < 0.1) galaxy clusters. 

Constraining the anisotropy of galaxy orbits from kine- 
matic data is a difficult task. The main problem arises 
from the so-called mass-anisotropy degeneracy which oc- 
curs in the Jea ns analysis of the velocity dispersion profile 
l|Merritt I [l987). In order to break this degeneracy a few 
solutions have been proposed. iBinnev fc Mamonl |l982l ) de- 
rived the velocity anisotropy profile by combining the Jeans 
standard analysis with an independent constraint on the 
mass profile. This approach was adopted by several authors 
who used the estimate of the mass profile from the kine- 
matics of elli pticals whose orbits app ear to be distributed 
isotropicall y (Biviano & Katgert 2004), lensi ng observations 
dNataraian fc Kneibl 1 19961 ) or X-ray data (jBenatov et al.l 
120061 ; iHwang fc Led I2OO8I ). Another method to solve the 
problem of the mass-anisotropy degeneracy relies on mod- 
elling t he degree of non-gaussia nity in the velocity distri- 
bution l|Merrifield fc Kentill990l ). A way t o account for this 
effect is to use two even v elocity moments l|Lokas fc Mamonl 
120031 ; lLokas et all |2006| ). Apart from this, in a few at- 
tempts a full model of the distribution function has been 
used, b ut in all cases a flat anisotropy profile was assume d 
(see e.g. Ivan der Marel et al . 2000; Ma hdavi fc Gellerjl2004l) . 
Nevertheless, this m ethod is still promising. As shown by 
IWoitak et ail (|2009h in a series of tests on mock kinematic 
data, the projected phase-space analysis wi th a suitable 
model of the distribution function proposed bv lWoitak et al.l 
(2008) allows to place robust constraints on the spatial vari- 
ation of the orbital anisotropy within the virial sphere as 
well as on the total mas s profile. In this pap er, we follow 
the method described in IWoitak et al.l l|2009l ) to constrain 
the orbital anisotropy at the cluster centres and around the 
virial radius. 

Most measurements of the orbital anisotropy in galaxy 
clusters have been obtained under the assumption of a flat 
profile. All res ults are consistent within errors with isotropic 



orbits (see e.g. van der Marel et al. 2000l; IRines et all 20031: 



Mahdavi fc Gellerjl2004l; lLokas et alj|2009: IWoitak fc Lokasl 



20071 ). Several studies were devoted to constraining the 
spatial variation of the orbital anisotropy. The analysis of 
galaxy kinematics of individual clusters points to rather 
isotropic orbits, although the anisotro py exhibits a signif- 
icant scatter between clusters (see e.g. iNataraian fc Kneibl 
1 19961 : iBenatov et"aT1l2006l ; IHwang fc Ledl2008l ). More robust 
constraints on the anisotropy profile were obtained in the 
analysis of the composite clusters resulting from combining 
a number of kinematic data sets into one. It was demon- 
strated that the orbits of early type galaxies tend to be 
isotropic within the virial sphere, whereas late type galax- 



ies are on predominantly radial orbits at t he virial radius 
l|Adami et al.lll998l ; lBiviano fc Katgertll2004l ) . This property 
is referred to as the orbital segregation and is commonly at- 
tributed to the orbital evolution from radially dominated 
orbits for spirals to isotropic ones for ellipticals. We note, 
however, that the orbital structure of both morphological 
types of galaxies can hardly be reconciled with the result of 
global isotropy. This problem is also addressed in this paper. 

The paper is organized as follows. In section 2, we de- 
scribe our cluster sample as well as kinematic data and the 
method of interloper removal. In section 3, we outline the 
details of our method of the projected phase-space analy- 
sis. The results of the analysis are presented in section 4. 
In the same section, we discuss the constraints on the mass- 
concentration relation and study the mass scaling relations. 
Section 5 is the devoted to comparison between our con- 
straints on the parameters of the mass profile and those 
found in the literature and based on other methods. In sec- 
tion 6, we present two tests of consistency between the data 
and some assumptions underlying our analysis. The sum- 
mary follows in section 7. In this work, we adopted a flat 
ACDM cosmology with density parameter f2 m = 0.3 and 
the Hubble constant Ho = 70 km s _1 Mpc - . 



2 DATA 

2.1 Cluster selection 

The primary criterion for cluster selection is the number 
of galaxies with spectroscopic redshifts per cluster. An opti- 
mum choice of a threshold for this number should result from 
the consideration of two contradictory effects. On one hand, 
the quality of statistical reasoning requires the number of 
redshifts to be as large as possible, but on the other we wish 
to build a sample with a statistically significant number of 
clusters which means that we also have to include clusters 
with a rather low number of spectroscopic redshifts. Trying 
to find a compromise we checked the number of available 
redshifts per cluster, as provided by NASA/IPAC Extra- 
galactic Database (NED). We found that the optimal choice 
for the purpose of our analysis is to consider clusters with 
at least 70 galaxies within the aperture of 2.5 Mpc and with 
velocities from the range ±4000 km s _1 in the rest frame of 
the cluster. The radius of the aperture equals to the virial 
radius of massive galaxy clusters and its role is to truncate 
all clusters to the approximate size of the virial sphere (in 
projection). The velocity range ±4000 km s _1 is a commonly 
used velocity cut-off which allows to perform a preliminary 
separation of cluster members from the interlopers (galaxies 
of background or foreground) . 

The second criterion for cluster selection concerns their 
dynamical state. The reason for considering this property is 
the fact that the basic assumption underlying our method of 
dynamical analysis is the state of equilibrium. It means that 
our sample should be possibly free of dynamically disturbed 
clusters, especially major mergers. In order to identify such 
clusters we studied X -ray images of c lusters obtained with 
ASCA llHorneijl200ll). XMM-Newto n (|Snowden et al.ll2008l ) 
and ROSAT (|Schuecker et aill200ll ) satellites. We rejected 
all clusters whose images exhibit signatures of merger activ- 
ity in the form of clearly visible distortions (e.g. A3667) or 
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the presence of massive subclusters (e.g. A3558, the most 
massive cluster in the Shapley Concentration). 

Another approach to assess the dynamical state of clus- 
ters was to study the general properties of galaxy kinemat- 
ics. This method was particularly important in the case of 
clusters for which X-ray data were not available. Following 
this approach, we excluded all clusters with clearly asym- 
metric or bimodal velocity distributions which are good sig- 
natures of a merger or at least an accidental alignment of 
some structures. Let us note that the analysis of galaxy kine- 
matics was carried out after the proper removal of interlop- 
ers (see subsection 2.4). 

As an independent test, we also checked for the 
presence of a cool core which ap pears to correlate 
with the dynamical state of c l usters |Burns et all 120081 ; 
ISanderson. Edge fe Smitbl 120091 ; iHudson et al l |2010| ) and 
may be regarded as an indicator of how relaxed the clus- 
ter is. We note, however, that the presence of a cool core 
was a secondary criterion of our selection. We used this cri- 
terion only in cases when the conclusions drawn from the 
kinematics or the X-ray image seemed to be ambiguous or 
contradictory. 

Finally, we also excluded clusters which appear to be 
accompanied by other clusters at comparable redshifts and 
relative positions on the sky smaller than 5 Mpc which is 
around twice the expected virial radius for a typical cluster. 
The majority of them are binary clusters at the initial state 
of merging, e.g. A399 and A401. The main reason for reject- 
ing these systems is the fact that their velocity diagrams are 
severely contaminated by galaxies from neighbouring clus- 
ters. Due to the proximity of these clusters both on the sky 
and in redshift space, velocity diagrams overlap and the se- 
lection of cluster members becomes impossible or at least 
highly ambiguous. 

2.2 Cluster sample 

In order to find all nearby (z < 0.1) galaxy clusters which 
satisfy selection criteria described above we searched NED 
and the recently published, but not yet included in NED, 
spectroscopic data from WINGS (The Wide field Nearby 
Galaxy cluster Survey) which is a photometric and spec- 
troscopic survey of a few tens of nearby galaxy clusters 
l|Cava et al.l 120091 ). The final sample resulting from our 
data mining consists of 41 galaxy clusters, o ut of which 
40 are rich clusters from the Abell catalogue (|Abelll 1 19581 ; 
lAbell. Corwin fc Olowfnll 19891 ) a nd one (MKW 4) is a poor 
clust er containing a cD galaxy (|Morgan. Kavser fc White] 
1 19751 ). Table [T] provides the list of all clusters as well as 
their basic observational properties, such as the heliocentric 
redshift z, the number of spectroscopic redshifts within the 
aperture of 2.5 Mpc and the velocity range ±4000 km s" 1 
around the mean JVtot and after the removal of interlopers 
(see subsection 2.4) A^mem- 

All clusters in Table [1] are relaxed objects in terms of 
merger activity. However, we cannot expect that the whole 
cluster sample is completely uniform in terms of the evolu- 
tionary stage. In order to check this, we searched for two 
indirect indicators of cluster evolution. The fi rst indicator 
is the pres ence of a coo l core. As shown by iBurns et al.l 
(|2008h and lHudson et al.l (|2010h . the survival of a cool core 
in nearby galaxy clusters is likely associated with the quies- 



Name 


z 


JVtot 


-^mcm 


B-M type 


cc 


A0085 


0.0551 


325 


317 


I 


yes a 


A0119 


0.0442 


215 


211 


II-III 


no a 


A0133 


0.0566 


75 


75 


I-II 


yes a 


A0262 


0.0163 


111 


135 


111 


yes a 


A0376 


0.0484 


85 


85 


1-11 


? 


A0496 


0.0329 


324 


313 


I 


yes a 


A0539 


0.0284 


141 


126 


111 


yes b 


A0576 


0.0389 


241 


214 


III 


yes a 


A0671 


0.0502 


129 


124 


II-III 


no d 


A0779 


0.0225 


151 


139 


I-II 


yes d 


A0954 


0.0932 


68 


66 


I-Il 


? 


A0957 


0.0436 


138 


124 


I-II 


no d 


A0978 


0.0544 


100 


90 


11 


? 


A1060 


0.0126 


369 


365 


111 


yes a 


A1139 


0.0398 


147 


143 


III 


? 


A1190 


0.0751 


137 


122 


11 


? 


A1314 


0.0335 


110 


127 


III 


no d 


A1650 


0.0838 


212 


204 


1-11 


yes a 


A1691 


0.0721 


118 


110 


11 


? 


A1767 


0.0703 


152 


147 


11 


no d 


A1773 


0.0765 


109 


103 


111 


? 


A1795 


0.0625 


177 


163 


I 


yes a 


A1809 


0.0791 


131 


123 


11 


no d 


A1983 


0.0436 


117 


100 


III 


yes d 


A2052 


0.0355 


150 


107 


I-II 


yes a 


A2063 


0.0349 


150 


115 


II 


yes a 


A2107 


0.0411 


88 


85 


I 


yes c 


A2142 


0.0909 


229 


226 


11 


yes a 


A2175 


0.0951 


93 


87 


11 


? 


A2415™ 


0.0581 


98 


96 


III 


yes d 


A2593 


0.0413 


230 


210 


11 


no d 


A2634 


0.0314 


223 


185 


11 


yes a 


A2670 


0.0762 


246 


238 


I-II 


yes c 


A2734 


0.0625 


174 


119 


III 


no b 


A3158 


0.0597 


148 


145 


I-II 


no a 


A3571 


0.0391 


180 


168 


I 


yes a 


A3581 


0.0230 


82 


74 


I 


yes a 


A3809™ 


0.0623 


105 


99 


III 


? 


A4059 


0.0475 


237 


192 


I 


yes a 


AS805 


0.0139 


159 


148 


1 


? 


MKW4 


0.0200 


132 


124 


1 


yes a 



Table 1. Observational properties of galaxy clusters selected for 
the analysis. The table includes the cluster name, heliocentric 
redshift (z), number of galaxy redshifts within the aperture 2.5 
Mpc and the velocity range ±4000km s" 1 (Ntot), number of red- 
shifts after interloper removal (-/V mem ), Bautz-Morgan morpholog- 
ical type (B- M type) and the com me nt on a cool core j ud ged on 
the ba sis of: ^Hudson et all told ) , ^Chen et all ||20071 ), IWhitd 
||2000|) . < iWmte, Jones fc Formanf Jl997l) . Superscript ui marks 
clusters from the WINGS survey. 



cent phase of recent history, when the cluster has not expe- 
rienced major mergers. Therefore the clusters which possess 
cool cores are expected to be more relaxed. The second indi- 
cator is the m orphological type. We u sed the Bautz-Morgan 
classification jBautz fc Morgan! [l970l ) which describes the 
degree to which a cluster is dominated by the BCG (bright 
central galaxy). Like for many other morphological classifi- 
cations of galaxy clusters, the Bautz-Morgan morphological 
type appears to correlate w ith the dynamical state of clus- 
ter evolution ( Sarazi nl ll988l ). The early type clusters (I,I-II) 
which are more dominated by BCGs are more dynamically 
relaxed than clusters of late types (II-III, III). 
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The last two columns of Table [T] include information on 
the presence of a cool core and the Bautz-Morgan morpho- 
logical type. Observational c onstraints on a cool core were 
obtained from the literature rtWhite. Jones k, Forman|[l997l ; 
IWhitell2000l ; IChen et al.ll2007l ; iHudson et al.ll20ld )~We note 
that the criteria used as the signatures of a cool core were 
not common to all authors and were mainly based on the 
inspection of such effects as a temperature drop at small 
radii, substantial mass deposition rate or sufficiently small 
cooling time in the cluster centres. Due to the fact that our 
cluster sample is optically selected, the X-ray data are not 
equally accessible for all objects. In particular, the informa- 
tion on the presence of a cool core was found for only 31 
clusters (22 cool core clusters and 9 non-cool core clusters). 
For other clusters X-ray data are very limited and therefore 
insufficient to study this property. 

The Bautz-Morgan morphological types of the clusters 
were obtained from NED. We note that this classification is 
not objective and may be treated as indicative only. Nev- 
ertheless, it allows us to divide our cluster sample into two 
classes, the early and late type clusters, which are thought to 
represent more and less advanced dynamical state of cluster 
evolution. A similar division of the clusters may be also car- 
ried out with respect to the presence of a cool core. In this 
approach cool core clusters are expected to be more dynam- 
ically evolved objects. A closer inspection of Table [1] shows 
that both classifications are clearly related to each other: 
the early type (I,I-II) clusters constitute 60 per cent of the 
cool core clusters and just 20 per cent of the non-cool core 
clusters. Let us note that this correlation, although not very 
tight, increases our confidence in using both classifications 
to assess the dynamical state of a cluster. 

2.3 Velocity diagrams 

The positions and redshifts of galaxies in the field of the clus- 
ters from Table Q] were obtained from NED. This database 
provides the most up-to-date compilation of many differ- 
ent kinds of extragalactic data, including spectroscopic red- 
shifts. The majority of spectroscopic data stored in NED 
comes from data releases of large scale redshift surveys, e.g. 
the S loan Digital Sky Survey (see lAdelman-McCarthv et al.l 



120081 . for the latest data release) or the 2dF Galaxy Red- 
shift Survey l|Colless et al.l [200ll ). as well as a number of 
spectroscopic surveys dedicated to observations of individ- 
ual g alaxy clusters, e.g. the NOAO Fundamental Plane Sur- 
vey (|Smith et alj|2004 ). We note that the SDSS, the 2dF 
survey and the NFPS are the three dominant sources of 
redshift data for our cluster sample. For two clusters, A2415 
and A380 9, redshifts and positions come from the WINGS 
survey l|Cava et al.ll2009T ). We found that in both cases the 
number of redshifts provided by the catalogue of this survey 
is considerably higher than that provided by NED. 

In order to calculate the projected clustercentric dis- 
tances R of all galaxies one needs to find cluster centres. 
From the observational point of view, the cluster centre may 
be defined as the position of the X-ray emission peak or the 
location of the BCG. For clusters that are not dynamically 
disturbed both centres coincide with each other and both 
point very well to the minimum of the gravitational poten- 
tial well. Since our cluster sample was constructed using op- 
tical data, we decided to assign cluster centres as the BCG 



positions. We found that 39 clusters from the sample have 
a well-confirmed central galaxy classified in NED as a cD 
galaxy (27 clusters) or a BCG (12 clusters). For two clus- 
ters, A2415 and A3809, central bright galaxies were found 
by visual inspection of optical images. Both clusters have 
distinct, bright and extended ellipticals lying in the close 
vicinity of the cluster centre given by the catalogue of Abell 
clusters. 

With cluster centres given by the positions of BCGs, 
we calculated the projected clustercentric distances of all 
galaxies. Except of a few clusters, the distances extend to 
J?max = 2.5 Mpc which is a fixed size of the aperture in 
our analysis. Three exceptional clusters are A376, A2415 
and A3809. The velocity diagram of the first of these is lim- 
ited to the aperture of around 1.8 Mpc which may be the 
effect of the incompleteness of the spectroscopic survey or 
may reflect extreme compactness of this cluster. Galaxy dis- 
tances of the other two clusters are limited by the aperture 
(aro und 1.5 Mpc) to which the WINGS survey is complete 
(see lCava et al.ll2009l ). For the analysis of these three clusters 
we assumed i? max = 1-5 Mpc. We note that an independent 
test for data completeness is presented in subsection 6.1. 



2.4 Interloper removal 

In order to separate cluster galaxies from the interlopers, i.e. 
galaxies of b ackground or foreground, we a pplied the method 
proposed bv lden Hartog fc Katgertl (1 19961 ) which appears to 
be o ne of the most effect i ve algorithms for interl oper removal 
fsee lWoitak et al.|[2007l ; IWoitak fc Lokas! 120071 ). In this ap- 
proach the interlopers are identified as galaxies with velocity 
along the line of sight exceeding a maximum velocity evalu- 
ated for a given projected radius R. Calculation of the max- 
imum velocity v ma x(R) is based on a toy model of cluster 
dynamics which assumes that galaxies follow circular orbits 
with velocity v c ir = \J GM(r)/r or fall towards the cluster 
centre with velocity V2v c i T . Evaluation of Vm ax (R) in such a 
model is simplified to finding the maximum of both velocity 
vectors projected onto the line of sight. We note that the 
interloper removal is an iterative procedure, i.e. the scheme 
is repeated several times until convergence is achieved. The 
mass M(r) in e ach iteration is app roximated by the virial 
mass estimator (|Heisler et al.lfl985l ) evaluated for all galax- 
ies within the aperture R = r. 



3 OVERVIEW OF THE METHOD 

In this section, we summarize our method of statistical in- 
ference of the mass and anisotropy profiles fro m galaxy kine- 
matic s, as developed and described in detail in lWoitak et al.l 
( 2009). The method relies on the Markov Chain Monte Carlo 
(MCMC) analysis of galaxy distribution in the projected 
phase-space, i.e. the space of the projected radii R and ve- 
locities along the line of sight vi os , in terms of a properly 
parametrized model of the distribution function. The pro- 
jected phase-space density fios(R,V\ os ) is related to the dis- 
tribution function f(E, L) through the following integral 



/losCR, wios) = 2tyR 



dz 



dv R dv,pf(E,L), 



(1) 



B>0 
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where z is the distance along the line sight from the cluster 
centre, vr and v$ are velocity components in cylindrical co- 
ordinates, E is the positively defined binding energy per unit 
mass, L is the absolute value of the specific angular momen- 
tum and z max is the distance of equality between vi os and 
the escape velocity for a fixed R. We note that the assump- 
tion underlying the functional form of both fi os (R, vi os ) and 
f(E, L) is spherical symmetry and the state of equilibrium. 

In our analysis we adopted the model o f the distribution 
function proposed bv lWoitak et al.l (|2008f ). The model was 
adjusted to the phase-space properties of dark matter par- 
ticles in cluster-size haloes formed in cosmological iV-body 
simulations. It assumes separability in the energy and the 
angular momentum, i.e. f(E,L) = /^(^/^(L), in which 
the angular momentum part is expressed by the following 
analytical ansatz 



tions 



-/3oo+/3 



-2/3 



(2) 



Po and Poo in ([2]) are two asymptotic values of the anisotropy 
parameter 



P{r) 



o*(r) 



(3) 



describing a ratio of the radial ay to tangential <jg veloc- 
ity dispersion. The profile of the anisotropy is a monotonic 
function chang ing from Po at the cluster centre to Poo for 
large radii fsee lWoitak et al1 l2008). The scale of transition 
between these two asymptotes is specified by parameter Lp. 

The energy part Je (E) is given by the integral equation 
defining the relation between the distribution function and 
the assumed density profile of a tracer p(r) 



p(r) 



-/3oo+/3 



-2/3(^3 



<Tv. 



(4) 



This equation can be simplified to the form of a one- 
dimensional int egral and solv e d num erically for /e(E) (see 
Appendix B in IWoitak et al.l I^POSt ) 1 ). In the following, we 
approximate the total density dis tribution by the universal 
NFW profile (jNavarro et al.Hl997t ). i.e. 



p(r/r s ) = 



47r(ln2- 1/2) (r/r s )(l + r/r s ) 2 rf ' 



(5) 



where r s is the scale radius and M s is the mass enclosed 
in a sphere of this radius. We also assume that galaxies 
trace dark matter. This assumption, commonly referred to 
as a hypothesis of constant mass-to-number density profile, 
was shown to be fairly held for relatively large clustercentri c 
radii (see e.g. lCarlberg et al ] |l997l ; iBiviano fc Girardill2003h . 
On the other hand, its robustness is more uncertain in 
the cluster centre: observations suggest that the profile 
of galaxy dis tribution at small radii ma y be as steep as 
dark matter (|Lin. Mohr fc Stanfordll2004r ) or exhibit a core 



l|Popesso et al. 2007 ). The problem of the consistency be- 
tween the assumption about mass-to-number ratio and the 
data collected for our analysis is addressed in section 6. 

Although M s and r s set the most natural scales for the 
phase-space units and both are used as the primary param- 
eters of the mass profile in our analysis, a common way 
to parametrize the universal mass profile is to use the virial 
mass M v and the concentration parameter c defined by equa- 



(4/3)7rr-3 



A c p c 

r v /r s 



(6) 



where r v is the virial radius, p c is the present critical density 
and A c is the virial overdensity. Fo r the cosmological mode l 
assumed in this paper A c = 102 |Brvan fc Norman] 1 19981 : 
lLokas fc Hoffmanbooir ). 

Constraints on the parameters of the mass and 
anisotropy profiles were determined by means of the analysis 
of the posterior prob a bility with the MCMC technique. Fol- 
lowing lWoitak et al.l (|2009l ) we used the likelihood function 
defined by 



L OC Y[ flas(Ri; 



cio 



\M S 



(7) 



where i is the reference number of a galaxy. The pro- 
jected phase-space density Q was evaluated numerically 
usi ng the algorithm of Gaussian quadrature, as outlined 
in IWoitak et alj (120091 ). We note that /i os preserves the 
normalization within the aperture 7i ma x fixed in subsec- 
tion 2.3. For the scale parameters we used the Jeffreys pri- 
ors, i.e. p(r s ) oc l/r s and p{M s ) oc 1/M S . In order to as- 
sign equal weights for all types of the orbital anisotropy we 
adopted P — s> (— l/2)ln(l — /3) = \n(a r /ae) reparametriza- 
tion and assumed flat priors for redefined parameters of 
the anisotropy profile. We fixed the value of Lq parame- 
ter at 0.45r s \J GM s /r s which corresponds to ~ lr s scale 
of the transition between /3o and Poo, as found for mas- 
sive dark matter h aloes from cosmological simulations (see 
IWoitak et al.ll2o"o"8h . In order to keep positivity of the phase- 
space density in the cluster centre we narrowed the range of 
the p rior for the central anisotropy to /3o < 1/2 (lAn fc Evansl 
2006). We also placed an upper limit on the outer anisotropy, 
Poo = 0.99, which solves the problem of an improper poste- 
rior distribution for models with Poo ~ 1- 

We determined the Markov chains using the Metropolis - 
Hastings algorithm (see lGregorvll2005l ; IWoitak et al.ll2009h . 
All free parameters of the algorithm were properly ad- 
justed to keep the acceptance rate (the probability of pa- 
rameter change between each two neighbouring models of 
the Markov chain) at the level of ~ 30 per cent which 
is the recommende d value for many-parameters models 
l|Gelman et alj|2004h . Our Markov chains consist of 2 x 10 
models. We checked that this length is sufficient to prop- 
erly explore our parameter space, also in a joint modelling 
described in subsection 4.5. The convergence of the MCMC 
algorithm has been tested by means of monitoring the profile 
of the posterior probability along the chains, the analysis of 
the autocorrelation functions calculated for all parameters 
in use and the comparison of the parameter dispersions eval- 
uated in different sections of the Markov chains. In this ap- 
proach, properly mixed Markov chain is characterized by the 
posterior probability fluctuating around a flat profile, fast 
decay of the autocorrelation functions and possibly small 
variations of the parameter dispersions. 
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Figure 1. Results of the MCMC analysis of the velocity diagram for the A1060 galaxy cluster. Black solid lines are the contours of 
the la and 2a credibility areas (panels below the diagonal) or the profiles of the marginal probability distribution (panels along the 
diagonal). Gray areas indicate the la credibility ranges of the parameters. 



4 RESULTS OF THE ANALYSIS 

We analyzed velocity diagrams of 41 nearby galaxy clusters 
in terms of the projected phase-space density, as described 
in section 3. Fig.[T]shows an example of the results obtained 
for the richest velocity diagram of the A1060 galaxy cluster. 
The contours show the la and 2a credibility constraints on 
the parameters of the mass profile, M s and r a , and the pro- 
file of the orbital anisotropy, j3o and /3oo . The panels along 
the diagonal and in the top right corner show the profiles 
of the marginal probability distribution and indicate the la 
credibility range for all parameters in use. We note that the 
la and 2a credibility ranges or contours are defined as those 
containing respectively 68 or 95 per cent of the correspond- 
ing marginal probability. 

To save space we do not show graphical representations 
of the results for other clusters. Instead, in Table [2] we list 
numerical constraints on the parameters in the form of MAP 
(maximum a posteriori) values and the la credibility ranges. 
Apart from the primary parameters of the mass profile used 
in our analysis, i.e. M s and r s , Table[2]also shows constraints 
on the virial mass M v and the concentration parameter c. 

Since our model of the distribution function is well- 
defined only within the virial sphere, we decided to replace 
the parameter /3oo with f3 Tv which is the anisotropy at the 



virial radius evaluated from the MAP value of the corre- 
sponding virial mass. The main motivation for this parame- 
ter replacement is to avoid the problem of the extrapolation 
of our model beyond the virial sphere. The new parame- 
ter of the orbital anisotropy was calc ulated using formula e 
for the velocity dispersions outlined in lWoitak et al.l (2009). 
The errors were evaluated by propagating the la limits for 
the asymptotic values of the anisotropy profile to the virial 
radius. The final constraints on the anisotropy parameters 
are summarized in Table [2] in terms of a commonly used 
anisotropy parameter, /3, as well as in terms of the ratio of 
the radial-to-tangential velocity dispersion, a r /a g . 



4.1 Scaling parameters of the mass profile 

Fig. [2] shows the constraints on the scale radius r s and the 
scale mass M a , as listed in Table [2] Both parameters ex- 
hibit a clear correlation which may be attributed mostly to 
the presence of a typical scale of the overdensity inside the 
sphere of the scale radius or, equivalently, to a typical value 
of the concentration parameter. In order to show this prop- 
erty more explicitly we plotted the line of constant c equal to 
7 which is the median value of the concentration parameter 
for the clusters (see Fig. [3]). This concentration parameter 
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Table 2. Constraints on the parameters of the mass and anisotropy profile of 41 galaxy clusters from the MCMC analysis of their 
velocity diagrams. The table provides the MAP values and the ranges containing 68 per cent of the corresponding marginal probability. 



corresponds to the mean density inside the sphere of scale 
radius equal to ~ 5600p c - 

Parameters of individual clusters are scattered around 
the line of the mean overdensity. This scatter is also shown 
in the form of the distribution of the concentration param- 
eter (see Fig. [3]). Interestingly, clusters with well-confirmed 
cool cores as well as clusters of early morphological type 
tend to populate the left part of the diagram in Fig. [2] This 
means that the clusters which are dynamically more evolved 
exhibit higher central overdensity (or higher concentration 
parameter). This property is also well-seen in Fig. [3l where 
we plotted the relative number of both types of clusters in 
subsequent bins of the concentration parameter. Cool core 
clusters (early type clusters) contribute 70 (60) per cent to 
the number of clusters with the concentration above the me- 



dian value and 40 (30) per cent to the number of clusters 
below the median value. Although the effect is not promi- 
nent, we can conclude that clusters with high concentrations 
are more likely to be dynamically more evolved than clusters 
with less concentrated mass distribution. 

The correlation between the state of dynamical evolu- 
tion of clusters and the steepness of the mass profile shown 
in Figs. [2] and [3] is in agreement with the results of recent 
cosmological simulations which point to two main factors 
responsible for this effect. First, more evolved clusters have 
not experienced severe mergers in their recent history and 
therefore tend to posses s dark matter haloe s with a steeper 
density profile (see e.g. iMaccio et al.l 120081 '). Second, quies- 
cent phase in the recent formation history of t hese clusters 
likely leads to the survival of cool cores (see e.g. iBurns et al.l 
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Figure 2. Parameters of the mass profile obtained from the 
MCMC analysis of 41 galaxy clusters, as listed in Table [2] Filled 
circles and crosses indicate early type (I,I-II type in Bautz- 
Morgan classification) clusters and cool core clusters. The solid 
line represents a family of density profiles with c = 7 which is 
the median value of the concentration parameter in the cluster 
sample. 



Figure 4. The mass-concentration relation from the analysis of 
41 galaxy clusters. The solid line shows the best fit power-law 
profile for the mean concentration parameter. Shading indicates 
logarithmic dispersion of the concentration distribution, In A, and 
the error of the normalization of the mean mass-concentration 
relation, c c (jvf v ). 
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Figure 3. The distribution of the concentration parameter from 
Table [2] (shaded area). Dashed and solid lines show the relative 
numbers (right axis) of the cool core and early type (I, I-II type 
in Bautz-Morgan classification) clusters, respectively. 



120081 : iHudson et alj|201(ih and triggers the accumulation of 
baryonic gas due to the process of radiative cooling (see e.g. 
iGnedin et alll2004h . Higher concentration of gas in cluster 
centres deepens the gravitational potential well and conse- 
que ntly makes the density profile of dark matter steeper (see 
e.g. iDolag et al.ll2009l ; iDuffv et alj|20ich . 



4.2 The mass-concentration relation 

The concentrati on parameter and the virial mass are weakl y 
anticorrelated (jNavarro et all 1 19971 ; iBullock et~ai] l200l] ). 
This relation reflects the mass assembly histories of dark 
matter haloes and may be used as an independent test of 
cosmological models. Here we use our constraints on the 
virial mass and the concentration parameter to test theoret- 
ical mass-concentration relation against observations in the 
mass regime of galaxy clusters. 



Fig. [4] shows our estimates of the virial mass M v and 
the concentration parameter c. According to the results of 
cosmological simulations the scatter of the concentration pa- 
rameter tnay_j2e_effiertively modelled by a logno rmal distri- 
bution (|Bullock et al.ll200ll ; lMacci6 et all 120081 ). The mean 
value of this distribution weakly varies as a power-law func- 
tion of the virial mass, whereas the logarithmic dispersion 
in the mass range of gal axy clusters is st atistically consis- 
tent with being constant (|Neto et al.ll2007| ). This means that 
the distribution of clusters on the mass-concentration dia- 
gram ma y be described by the following probability density 
function |Comerford fc Nataraiar]|2007l ) 



/( C )dC : 



1 



27rc In A 



exp 



(In c — In fxf 



2 (In A) 2 



dc, 



(8) 



where fi = c(Mq){M v / 'Mo) a . In order to minimize the cor- 
relation between the slope and the normalization, the ref- 
erence mass Mo should be approximately the median value 
of the virial mass in the sample. In our case we assumed 
M = 5 x 10 14 M Q . 

We carried out the likelihood analysis of the data shown 
in Fig. 13] The likelihood function was the product of the 
probability density function (|8]) evaluated at the mass pa- 
rameters of all clusters. The analysis was performed with the 
MCMC technique with three free parameters: the slope a, 
the normalization c(5 x 1O 14 M0) and the logarithmic disper- 
sion In A. We note that the latter parameter includes both 
the internal scatter of the theoretical mass-concentration re- 
lation and the inaccuracy of the measurement of the concen- 
tration parameter. 

Fig. [S] shows the constraints on the slope and the nor- 
malization of the mass-concentration relation. Three con- 
tours correspond to the la, 2a and 3a credibility regions. 
Empty symbols with error bars indicate the constraints on 
the slope and the normalizati on of the M v — c relation ob- 
tained by iBuote et all (l2007h from the a nalysi s of X-ray 
clusters and by IComerford fc Nataraianl ((2007) from the 
analysis of the mass parameters of galaxy clusters compiled 
from the literature. In both cases, the original results were 
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Figure 5. Constraints on the parameters of the power-law fit 
to the mass-concentration data, the slope a and the normal- 
ization at the mass 5 X 1O 14 M0. The empty square and circle 
with error bars show the constraints o n the parameters o f the 
mass-concentration relat i on ob tained bv lBuote et al.l J2007T) and 
IComerford fc Nataraianl 1120071 ). respectively. Filled symbols in- 
dicate the calibrations of the mass-co ncentration relation from 
cosmological TV-body simulations (see IMaccio et al.l l200Sh with 
the parameters from three data releases of the WMAP satellite 
labelled by WMAPI, WMAP3 and WMAP 5. 



properly adapted to the overdensity parameter assumed in 
our work as well as to the reference mass Mo. Our re- 
sult agrees within erro rs with the constraints obtaine d by 
iBuote et ail l|2007h and IComerford fc Nataraianl i|2007f ). al- 
though clearly a lower value of the normalization is pre- 
ferred, c(5 x 10 14 M Q ) = 6.9±° - 6 7 . 

In order to compare our results with theoretical pre- 
dictions we plotted in Fig. [5] the calibrations of the mass- 
concentration relation for relaxed dark matter haloes formed 
in c osmological TV-body simulations (see IMaccio et al.l 
2008), where an appropriate change of the virial overdensity 
parameter was taken into account. We note that the refer- 
ence simulations are purely dark matter simulations which 
do not take into account the steepe ning of dark matter pro- 
file induced by baryonic gas (see iGnedin et ail I2004T ), To 
account for this effect we increased all normalizations by 
10 per cent which is the value close to the upper limit of 
the bias in the concentration parameter foun d in hydro- 
dynamical simulations of galaxy clusters (see iDolag et al.l 
120091 ; iDuffv et al.|[2O10l ). Each calibration (filled points in 
Fig. refers to one of three cosmological models given by 
the parameters from three data releases of the Wilkinson 
Microwave Anisotro py Probe ( WMAP) satellite, W MAPI 
l|Spergel et ail 120031). WMAP3 (|Spergel et al.l |2007| ). and 
WMAP5 (|Komatsu et al.l 120091 ). The main difference be- 
tween these three sets of cosmological parameters lies in 
the value of as, the linear amplitude of density pertur- 
bations at 8 /i _1 Mpc scale at 2 = 0, which ranges from 
0.75 for WMAP3, through 0.8 for WMAP5, to 0.9 for 
WMAPI. Keeping in mind that the normalization of the 
M v — c relation is mostly sensitive to as, the three val- 
ues of as correspond to three normalizations labelled by 
WMAPI, WMAP3 and WMAP5 in Fig-.g] From Fig.[S]one 
can see that our constraints on the normalization of the 
mass-concentration relation are fully consistent with as of 



WMAPI cosmology and agrees only at 2a and ~ la level 
with that of WMAP3 and WMAP5, respectively. We con- 
firm recent findings that the concentration parameters for 
galaxy clusters tend to be higher than those predicted by 
cos mological model w i th WMAP3 or W MAP5 parameters 
(see lBuote et al.ll2007l ; IDuffv et al.ll2008l ). 

In order to place constraints on as using the lim- 
its for the normalization c(Mo) we adopted a semianalyt- 
ical model for the m ass-concentration relation proposed by 
iBullock et al.l (|200ll ) . According to this model the mean con- 
centration parameter at redshift z — is inversely propor- 
tional to the scale factor a c of halo collapse given by the 
following equation 

where D(a) is the linear growth rate, a is the present lin- 
ear rms density fluctuation at the mass scale FM V and F 
is a free parameter of this model. Assuming that the only 
factor determining the normalization c(Mo) at a fixed virial 
mass Mo is the value of as we find that erg oc 1/D(a c ) and 
c(Mo) oc l/a c . We used both proportionalities to convert 
c(Mo) int o as- In this convers ion, we assumed the calibra- 
tion from IMaccio et al. I l|2008l) for as = 0.9 and corrected 
for 10 per cent bias upwards, i.e. c(5 x 10 14 Mq) = 6.8 (see 
Fig. [5j, where a c was calculated using equation <(9j with 
F = 0.001 which is the r ecommended value for the mass 
regim e of galaxy clusters (|Bullock et al.l l200ll ; IBuote et al.l 
l2007h . As a result we expressed the constraints on the mean 
mass-concentration relation in terms of as parameter ob- 
taining as = 0.91± o ;^. We note that the difference between 
MAP value of as and its value from WMAP5 data release 
corresponds to the change of c(Mo) by —14 per cent. We 
cannot exclude that this difference may arise from system- 
atic errors in the measurement of the concentration param- 
eter, observational selection biases or an insufficient value of 
the bias correcting the calibration of the mass-concentration 
relation for the presence of baryons. 

Our constraint on the slope of the M v — c relation is 
rather weak, i.e. a — —0.07 ± 0.15. It is consistent with a 
flat profile as well as with the slope obtained in cosmological 
simulations (a ~ —0.1). The best-fit value of the logarith- 
mic dispersion parameter In A is 0.55. The width of the light 
gray region in Fig. U indicates this value, whereas the black 
solid line and the dark gray area show the best-fit profile for 
the mean concentration /i(M v ) and the error of the normal- 
ization. Taking into account that the statistical scatter in 
the measurement of the concentratio n parameter is around 
a In A w 0.44 (see lWoitak et al.ll2009l ), we conclude that the 
logarithmic dispersion of an intrinsic scatter in the theoret- 
ical M v — c relation is about 0.35 which is in fair agreement 
with the results of cosmological simulatio ns (|Bullock et al.l 
120011 ; iNeto et al.ll2007l ; IMaccio et aill2008l ). 



4.3 Mass scaling relations 

The virial mass of galaxy clusters scales with a number of 
different observables. Here we consider scaling relations with 
the velocity dispersion, the X-ray temperature and the X-ray 
luminosity (see Fig.|SJ). The virial mass, as listed in Table [2J 
comes from our analysis of 41 velocity diagrams of nearby 
galaxy clusters. Velocity dispersions were calculated using 
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relation a M [10 14 M o ] 



M v 


= M 


a los /(750 km s" 1 ) 


3.17 ±0.13 


5.99 ± 0.17 


M v 


= M 


L x /(10 44 erg s" 1 ) 


0.49 ± 0.05 


5.53 ± 0.33 


M v 


= M 


-| a. 

T x /(5 keV) 


1.54 ±0.12 


7.85 ± 0.55 



Table 3. Best-fit parameters of the power-law approximation for 
the three different mass scaling relations. 



all galaxies lying within the virial radii of the clusters. The 
errors were estimated by bootstrapping from the sample. 
We used the same kinematic data as in the proper analysis 
of velocity di agrams. The temperatures were obtained by 
iHornerl (|200lh from the spectral analysis of data acquired 
by the ASCA satellite (23 clusters). The luminosities were 
taken from the catalogues of X-ray clusters observed by the 
ROSAT satellite in the energy band 0.1-2.4 keV: the North- 
ern RQSAT All-Sky (NORAS) galaxy cluster survey (see 
IBohringer et all l2000r i . the ROSAT-ESO Flux Limited X- 
ray ( REFLEX) galaxy cluster survey (see IBohringer et al.l 
2004|) and a comple mentary source of the data provided by 



Ledlow et ail ((2003). All luminosities were corrected for the 



assumed cosmological parameters. 

All scaling relations seen in Fig. |S]are well approximated 
by a power-law function. Table [3] lists the best-fit slope and 
normalization for each of them. Parameters of the scaling re- 
lations were obtained by fitting a line to the data in logarith- 
mic scales. We used symmetrized logarithmic errors given by 
[x<2 — xi)/(2xq), where xi and X2 are the boundaries of the 
lcr credibility range of x variable, and xo is the most prob- 
able value. The best-fit slope and normalization wer e cal- 
culated using the bisector method (|lsobe et al.| [l990). The 
pivot points of <7i os , Lx and T x given in Table |3] were taken 
to be approximately the median values of the corresponding 
variables. This minimized the correlation between the slope 
and normalization. The errors were evaluated by bootstrap 
resampling. 

The best-fit scaling relations were plotted in Fig. [6] 
Solid lines are the best power-law fits and the shaded ar- 
eas indicate the errors in both the slope and normalization. 
Both the temperature and the velocity dispersion appear 
to correlate very well with the virial mass. We note, how- 
ever, that the a tight correlation with the velocity disper- 
sion occurs mainly due to the fact that the virial masses 
were derived from the same kinematic data that were used 
to evaluate velocity dispersions. Two clusters, A119 and 
A2107, lie below the mean mass-temperature relation. This 
discrepancy is probably caused by the fact that these clus- 
ters are less relaxed compared to the other ones in our sam- 
ple: A 119 exhibits distinct substructering |Schuecker et al.l 
l200lf ), whereas A2107 has a high-velocity (~ 200 km s" 1 ) 
cD galaxy that may suggest some merging activity. The cor- 
relation between X-ray luminosities and the virial masses is 
weaker and confirms a well-known fact that the X-ray lu- 
minosity is not a robust mass proxy. Nevertheless, we note 
that our constraint on the slope of the M v —Lx relation is in 
fair agreement with other measurements avail able in the lit- 
erature ranging between 0.5 and 0.7 (see e.g. iPopesso et all 
l2005l ; lLopesetal .1120091 ). 
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Figure 6. Scaling relations between the virial mass and different 
mass proxies: the velocity dispersion (cri oa ), the X-ray luminosity 
(L x ) and the X-ray temperature (T x ). Solid lines show the best 
power-law fits and gray shaded areas represent the errors in both 
the slope and normalization. 



The slopes of the mass scaling relations with the X- 
ray temperature as well as with the velocity dispersion are 
fully consistent with the predictio ns of the virial theorem 
which states that M v tx a 3 tx T x /2 (|Brvan fc Normanlll99ct ). 
The slope of the mass-temperature relation is in excellent 
agreement with the most recent est imates which yield va l- 
ues between 1.5 and 1.6 (see e.g. IVikhlinin et al l I2006T I. 
The normalization of the same relation, when converted 
to the overdensity A c = 500, is (4.7 ± 0.3) x 10 14 M Q and 
agrees with that derived from high-resolution hydrodynami- 
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Figure 7. The orbital anisotropy of galaxies in 41 galaxy clusters. 
Empty and filled symbols refer to the clusters with the concentra- 
tion parameter above or below the median concentration in the 
sample (c = 7). Numbers indicate fractions of clusters in domains 
of the diagram defined by the lines of isotropic velocity distribu- 
tion (/3o = or /3r v = 0) and the family of flat anisotropy profiles 

03o = Av). 



cal simulations (see e.g. iBorgani et al.l l2004:) . It is consistent 
with the normalizat ion of X-ray derived mass-temperature 
relation (see e.g. Arnaud. Pointecouteau fc Prattj 1 2005; 
IVikhlinin et al.ll2006l ) and is lower by ~ 20 per cent than 
mass-temperature nor malization determined from weak 
lensing masses (see e.g. lPedersen fc Dahldl2007h . 

Our analysis provides an independent test of the scaling 
relations with X-ray observables on the scale of the virial 
overdensity that is usually not used for X-ray clusters. Since 
X-ray flux is detected in the very central part of clusters, 
a typical scale of the mass probed in X-ray observations is 
between M200 and Msoo- This mass scale is also commonly 
used to study X-ray scaling relations. Our analysis shows 
that the mass scaling relations both with Tx and Lx may 
be extended to the virial mass with the mean overdensity 
~ 100p c . 



4.4 The anisotropy of galaxy orbits: 
clusters 



individual 



Fig. [7] shows the constraints on the orbital anisotropy of 
galaxies in the cluster centre ( /3q) a nd at the virial sphere 
(/3 rv ). Although IWoitak et all (|2009h showed that the esti- 
mate of the anisotropy in the outer part of the cluster ap- 
pears to be much less accurate than in the central part, we 
think that some interesting conclusions may be drawn from 
studying the overall distribution of the results obtained for 
individual clusters. The estimates of the central anisotropy 
for all clusters are statistically consistent with /3o = 0. 
Galaxy orbits at the virial sphere of the clusters are no- 
ticeably radially biased with a typical anisotropy of around 
(cr r /ag) rv — 1.7. We also note that there is a clear tendency 
of the anisotropy profile to increase with radius: over 70 per 
cent of the results lie above the line of flat anisotropy profiles 
indicated by the diagonal. On the other hand, some clusters 
exhibit flat j3 profile with an isotropic velocity dispersion 
tensor. 



Clusters with different degrees of dark matter concen- 
tration tend to populate different regions of the parame- 
ter space shown in Fig. In particular, low concentration 
clusters (c < 7) exhibit more radially biased galaxy or- 
bits at the virial sphere than clusters with higher value of 
the concentration (c > 7). Furthermore, the former clus- 
ters seem to constitute a very homogeneous class in terms 
of the anisotropy profile. The orbital anisotropy of galax- 
ies in these clusters seems to grow with radius between two 
universal asymptotes: a r /ag ~ 1 in the cluster centre and 
ay/oa ~ 1.8 at the virial radius. In the case of the high 
concentration clusters, the anisotropy profiles do not show 
any signatures of such universality. The anisotropy at the 
virial radius in these clusters ranges from (a r /cre)r v ~ 1 to 
(cr r /ag) rv ~ 2. Nevertheless, it is interesting to note that 
85 per cent of the clusters for which constraints on the 
anisotropy parameters are statistically consistent with the 
fully isotropic model (/?o = Av = 0) are high concentration 
clusters. 

Keeping in mind that the degree of mass concentration 
in galaxy clusters is related to the state of dynamical evo- 
lution, we suggest that the dichotomy between the low and 
high concentration clusters observed in Fig. [7] may be likely 
the result of the cluster evolution. In this scenario galaxy 
orbits of less evolved clusters (low concentration parameter) 
would be dominated by radial motions in their outer parts 
(filled symbols in Fig. [7)l , whereas for more evolved clusters 
galaxy orbits tend to be more isotropic both in the cluster 
centre as well as in the vicinity of the virial sphere. 

At the early stage of cluster evolution one can expect 
an excess of radial orbits around the virial sphere. This 
excess is mostly the effect of radial infall of galaxies onto 
the clusters. Due to the relaxation process the initial orbits 
of infalling galaxies are expected to become more isotropic 
with time so that the most evolved clusters are character- 
ized by more isotropic distribution of galaxy orbits. A simi- 
lar mechanism of the is otropization of galaxy orb its in clus- 
ters was pointed out bv lBiviano fc Katgertl (|2004j ). Studying 
the orbital structur e in ENACS (ESQ Nearby Abell Clus- 
ter Survey) clusters, iBiviano fc Katgertl (|2004h showed that 
early type galaxies (ellipticals) are characterized by more 
isotropic orbits than spirals for which a r /ae at the virial 
sphere is around 1.7. Since the former galaxy population is 
more evolved and relaxed than the latter, this orbital segre- 
gation implies that a mechanism of the orbital isotropization 
must operate during cluster evolution. 



4.5 The anisotropy of galaxy orbits: a joint 
analysis 

In order to estimate the typical anisotropy profile of galaxy 
orbits in the clusters of our sample we carried out a 
joint MCMC analy sis of all velocity diagrams. Following 
IWoitak et ail (|2009h we defined the likelihood function as 

f(M s ,i, . • . , M s 

00 ) — 

n iVmom(j) 

II II fl*{Rj*,^jA{M.j,r.j,fa,P oa }), (10) 

j=l i=l 

where j and i are respectively the reference numbers of a 
clusters and a galaxy of the j-th velocity diagram. The two 
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Figure 8. Constraints on the parameters of the anisotropy pro- 
file of galaxy orbits from the joint analysis of 41 nearby galaxy 
clusters. The gray shading shows the lcr and 2a confidence re- 
gions of the marginal probability distribution and the cross indi- 
cates the MAP value. The black solid line marks the boundary 
of the parameter space excluded from the analysis by the prior 

{u r /crg)oc < 10. 

parameters of the anisotropy profile, /3o and poo, are com- 
mon to all clusters and define the anisotropy profile in clus- 
ter population, whereas parameters of the mass profile fix 
phase-space units for individual clusters. We used the con- 
straints on the parameters of the mass profile, M s and r s , 
in the form of Gaussian prior, where the means and the dis- 
persions were given by MAP values and symmetrized errors 
corresponding to the la credibility ranges. 

Fig. [8] shows the result of a joint analysis in the form of 
the la and 2a credibility contours of the marginal prob- 
ability distribution. In order to avoid extrapolation be- 
yond the virial sphere we converted /3oo parameter into the 
anisotropy at 7r s which is a typical scale of the virial ra- 
dius in the cluster sample. The black solid line indicates the 
boundary of the parameter space excluded from the anal- 
ysis by the prior (a r /ag)oo < 10. The joint analysis con- 
firms that galaxy orbits are typically isotropic in the cluster 
centres, (a r /ag)o = 0.97 ± 0.04, and radially anisotropic 
at the virial sphere, (oy / 'ae)i r ^ — 1.75lg'ig. The degree of 
the anisotropy at the virial sphere agrees with the value 
of the orbital anisotropy determined for late type galax- 
ies which dominate in t he outer parts of galaxy clusters 
jBiviano fc Katgertll2004). We do not confirm the result ob- 
tained by Ivan der Marel et~ai1 (|2000h who concluded that 
velocity distribution in galaxy clusters is consistent with 
isotropic orb its. The reason of this disc repancy may lie in 
that fact that Ivan der Marel et al.l (|2000T ) assumed a flat pro- 
file of the anisotropy. Taking into account that the central 
anisotropy is much better c onstrained than the anisotropy in 
the outer part of a cluster (|Woitak et al.|[2009T ). one may ex- 
pect that the estimate of the global anisotropy is dominated 
by its central value. Consequently, the final estimate of the 
anisotropy is closer to its central rather than the outermost 
value. 

The size of the credibility regions in Fig. [8] does not re- 
flect the scatter of the anisotropy profiles in the cluster sam- 
ple. It rather describes the accuracy of parameter estimation 
under the assumption that the cluster sample is homoge- 



Figure 9. Constraints on the parameters of the anisotropy pro- 
file of galaxy orbits from the joint analysis of low concentration 
clusters (with c < 7 and the median virial radius r v ~ 4r s ) and 
high concentration clusters (with c > 7 and the median virial 
radius r v ~ llr s . The gray shading shows the la and 2a con- 
fidence regions of the marginal probability distribution and the 
crosses indicate the MAP values. The black solid line marks the 
boundary of the parameter space excluded from the analysis by 
the prior (a r /ae)oc < 10. 

neous in terms of the anisotropy profile. In order to illustrate 
the degree of the internal scatter of the anisotropy parame- 
ters within our cluster sample we repeated the joint analysis 
for two separate cluster subsamples, clusters with the con- 
centration parameters above (c > 7) and below (c < 7) 
the median (see Fig. [9j . The choice of these subsamples is 
motivated by the observed segregation of low and high con- 
centration clusters on the plane of the anisotropy param- 
eters (see Fig. [7J). Our result confirms the conclusion from 
the previous subsection that low concentration clusters have 
galaxies on clearly radially biased orbits at the virial sphere, 
whereas galaxy orbits for high concentration clusters tend 
to be isotropic at all radii within the virial sphere. 



5 COMPARISON WITH OTHER METHODS 

For a number of galaxy clusters parameters of the mass pro- 
file have already been determined with other methods. In 
this section, we use these estimates to compare with the re- 
sults obtained with our method. This comparison allows us 
to check the credibility of our approach with respect to other 
methods of mass inference for galaxy clusters. 

Fig. [10] shows the comparison between the virial mass 
estimated in our analysis and with other methods avail- 
able for nearby galaxy clusters. The bottom panel in- 
cludes the mass estimates obtained from the X-ray obser- 
vations. It comprises both the r esults of a detailed analy- 
sis of the temperatu re profiles (|Pointecouteau et al.l 120051 ; 
IVikhlinin et~al1l2006t ) and the mass estimates obtained un- 
der the assumption that the intracluster gas is isothermal 
l|Reiprich fc Bohringerl [2002) or is linked to the tempera- 
ture via a polytropic equation of state |Sanderson et al] 
2003). The top panel presents the mass estimates from 
galaxy kinematics including the Jeans an alysis with the ve- 
locity dispersion <ti os and the kurtosis ki os l|Lokas et alfcOOol ; 
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Figure 10. Comparison between the virial mass obtained with 
the method of the projected phase-space analysis (M v ) and other 
methods based on the X-ray observations (M V ( X ) > bottom panel), 
the Jeans analysis of velocity moments and the caustic tech- 
nique (M v ( ve i), top panel). The sources of data for subsequent 
point types are t h e following: filled circl es in the top panel: 
lLokas et al (2006); Woitak & Lokas (2007 ); empty ci rcles in the 
top panel: iRines et al. (2003); Rincs & Diafcrio (2006); filled cir- 
cles in the bottom pan e l: iReiprich fe Bohringcr (20 021): empty 
squares: IVikhlinin et all j2006h; em pty tr iangle: Sand erson et all 
l|2003l l: cross: IPointecouteau et all (20051 ). The black solid line 
represents equality between the mass estimates. 



Woitak fe Loka1l2007l) and t he caustic method (|Rines et all 



Rines fc Diaferid l200fj| ). We note that all mass esti- 



mates were adopted to the virial overdensity assumed in 
this work. In the conversion between masses corresponding 
to different overdensities, we used the same mass profile as 
assumed by the authors of a given data analysis. Apart from 
the analysis carr i ed out bv lReiprich fc Bohringerl l|2002r ) and 
ISanderson et all (|2003h who used the mass profile derived 
from the equation of hydrostatic equilibrium with a fixed 
parametrization of the gas density and the temperature pro- 
file, the mass profile was parametrized by the NFW formula 
©■ 

Our mass determinations agree within errors with the 
estimates based on X-ray analysis. There are two clusters, 
A119 and A262, for which the difference between mass esti- 
mates exceeds significantly the la level. The reason of this 
discrepancy seems to be clear for A119, s ince this is one 
of th e least relaxed clusters of our sample l|Schuecker et all 
2001). On the other hand, an analogous discrepancy for 




Figure 11. Comparison between the concentration parameter 
obtained with the method of the projected phase-space anal- 
ysis (c) and from the X-ray observations (cx). The sources of 
data for subseq uent point type s are the following : filled squares: 
Xu et al. I ll200lh ; empt y squares: iGastaldello et all (1200711 : crosses: 
Vikhlinin et all (120061): circles: lEttori. De Grandi fc Molendil 
(120021) ; triangle: IPointecouteau et all (120051) . The black solid line 
represents equality between both estimates of the concentration 
parameter. The dashed lines connect data for the same clusters. 



A262 remains unclear, since there is no observational evi- 
dence which would indicate that this cluster is not in equi- 
librium. We checked that both the X-ray mass and the mass 
inferred from galaxy kinematics, as displayed in Fig. 1101 are 
consistent with t he corresponding result obtained by inde - 
pendent authors (|Lokas etalj|2006| ; IGastaldello et alll2007M . 
This suggests that this difference in probably a consequence 
of a generic tension between kinematic and X-ray data. 

The virial masses obtained in our work agree very well 
with those determined in the analysis of velocity moments. 
This fact is fully understandable, since our approach to data 
modelling is an extension of the formalism based on velocity 
moments. We note that the la credibility ranges obtained 
with our method are noticeably narrower than those result- 
ing from the Jeans analysis and are comparable to the errors 
of the X-ray masses. 

The caustic method provides mass estimates which 
marginally agree with the virial masses determined with 
our approach. The caustic mass correlates with the virial 
mass from our analysis, but is considerably scattered around 
the line of equal mass estimates and tends to be underes- 
timated with respect to our mass determination. For three 
cases, A85, A2142 and A2670, this effect is most prominent: 
the caustic virial mass is ~ 3 smaller than its counterpart 
from our analysis. We checked that our measurements of 
the virial masses for these clusters agree with the estimates 
from the X-ray observations yielding ~ 1.1 x 10 15 Ma for 
A85 (|Reiprich fc Bohringerl l2002i ; iDurret. Lima fc Forman 
l200rj) . ~ 2.1 x 10 15 M Q for A2142 (IReiprich fc Bohringei 
7 x 10 14 Af Q for A2670 (Sandcrs onet alj|2003h . 



2002T) and . 

This confirms that the caustic masses of these three clusters 
are likely underestimated. 

Fig. [TT] shows the comparison between the concentra- 
tion parameters obtained with our method and fr om the X- 
ray o b servations of the corresponding clust e rs (see Xu et all 
20011 ; lEttori. De Grandi fc Molendil |2002| ; IVikhlinin et"afl 
2009 : IGastaldello et al.ll2007l ). We note that some clusters 
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(oy/°e)o 

(tT r /o'e)oa 



1 Mpc 

0.001 ±0.140 
-0.030 ±0.171 
0.038 ±0.173 
0.025 ±0.105 
-0.035 ±0.101 
-0.112 ±0.588 



1.7 Mpc 

0.015 ± 0.056 

0.012 ± 0.094 
-0.010 ±0.101 

0.008 ± 0.042 
-0.008 ± 0.038 

0.001 ± 0.253 



0.017 ±0.049 
0.020 ± 0.077 
-0.019 ±0.082 
0.003 ± 0.033 
-0.003 ± 0.026 
-0.007 ±0.219 



Table 4. Statistics of the relative residuals of the model param- 
eters from the analysis of 41 velocity diagrams with three differ- 
ent truncation radii, i? ma x. Each row lists the means and disper- 
sions in the cluster sample for the relative residuals defined by 
log(Jf ij max /^2.5Mpc), where ^i{ max is the MAP value of the pa- 
rameter obtained in the analysis of velocity diagrams truncated 

at /?-max - 




1.0 

g o.o 

-0.5 
-1.0 



(A85, A119, A262, A496, A1795, A3571 and MKW4) have 
multiple representations on this diagram corresponding to 
the estimates obtained by different authors. Since for some 
clusters these estimates do not agree within the errors, we 
decided to plot all results. The points referring to the same 
cluster were connected with a horizontal dashed line. As for 
the comparison of the virial masses, the concentration pa- 
rameters were properly converted to the virial overdensity 
assumed in this work. 

The concentration parameter determined with our 
method correlates with that inferred from the X-ray data, 
although the scatter is significantly larger than for the virial 
mass. Assuming a power-law relation between both esti- 
mates, c = CY(cx/7) a , we obtained the following constraints 
on the slope and the normalization: a = 1.04 ± 0.4 and 
C7 = 8.2 ± 1.0, where the fit was performed in the same 
manner as described in subsection 3.3. Although the data 
sample is not statistically uniform so we cannot draw any 
firm conclusion, we point out that X-ray derived concentra- 
tions tend to be smaller by ~ 15 per cent than the concen- 
trations inferred from galaxy kinematics. 



6 SELF-CONSISTENCY TESTS 
6.1 Completeness 

An important factor that may influence the quality of our 
analysis is the completeness of spectroscopic data. In or- 
der to check this effect we carried out the following test. If 
the data are incomplete one may expect that the parame- 
ters estimated from velocity diagrams truncated at different 
projected radii will be different. To see if this is the case 
we repeated the analysis for all clusters with three different 
values of the maximum projected radius: _R max = 1 Mpc, 
-Rmax = 1-7 Mpc and i? max = r v , where r v is the virial ra- 
dius taken from the first analysis with _R max = 2.5 Mpc. 
Table |4] shows the comparison between parameter estimates 
obtained for different sizes of the aperture. Each row con- 
tains the logarithmic mean and dispersion of the relative 
residuals for the MAP values, where the reference values of 
the parameters come from the analysis with i? max = 2.5 
Mpc (see Table 0). 

We see that there is no trend in the mean value of 
the residuals. One may notice a negative bias in the outer 
asymptotic value of the anisotropy profile for i? max = 1 
Mpc. This is, however, a consequence of the fact that using 



Figure 12. Surface number density of galaxies in 41 galaxy clus- 
ters combined into one. The black solid line is the projected NFW 
profile. The vertical line indicates a typical scale of the virial ra- 
dius. 

a smaller aperture we probe the inner part of the clusters 
for which the orbits are more isotropic. We also note that 
the logarithmic dispersions which describe the scatter of the 
relative residuals are smaller or comparable to the relative 
errors associated with the la credibility ranges listed in Ta- 
ble [2] This property together with the stability of the mean 
values means that our results are not sensitive to the change 
of the truncation radius of velocity diagrams. This excludes 
the possibility that the data may be subject to a significant 
incompleteness effect. 

6.2 Mass-to- number density 

As a second test of self-consistency between the assumptions 
of our analysis and the data we show in Fig. [12] the sur- 
face number density of galaxies in all clusters combined into 
one. The surface density was evaluated in radial bins equally 
spaced in logarithmic scale, where the universal scaling of 
radii by r s was adopted (N s is the number of galaxies within 
the R < r s aperture). The points above 3r a were corrected 
for the incompleteness of velocity diagrams which occurs 
due to the fact that the value of i? ma x/»"s is not universal 
in the cluster sample. As a line of reference we plotted the 
projected NFW profile (black solid line). Apart from radii 
above ~ 3r s , where some deviations between the model and 
the data occur mostly due to the problem of incomplete- 
ness, galaxy density is well approximated by the NFW pro- 
file. Our data for small radii do not reveal a signature of a 
core i n galaxy distribution, as advocated bv lPopesso et all 
(2007). Although our analysis concerns the projected den- 
sities, we may conclude from Fig. [T5] that the assumption 
of our phase-space density model about a constant mass-to- 
number density is also satisfied. This confirms the observa- 
tional fa ct that has already been discussed in the liter ature 
(see e.g. ICarlberg et"ai1ll997l ; iBiviano fc Girardill2003r i. 



7 SUMMARY AND CONCLUSIONS 

In this paper we presented the analysis of 41 nearby (z < 
0.1) relaxed galaxy clusters. Using the method based on 
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the modelling of velocity diagrams in terms of the pro- 
jected phase-space density, as described in section 3, we con- 
strained parameters of the total mass and anisotropy profiles 
in these clusters. Parameter estimation was carried out in 
the framework of Bayesian data analysis with the use of the 
MCMC technique. 

We found that the concentration parameter is weakly 
correlated with the morphological type of the clusters as well 
as with the presence of a cool core. We interpret this fact as 
a signature of cluster evolution from the state with a flatter 
density profile (low concentration parameter) to one with 
a steeper density profile (higher concentration parameter). 
Our explanation is consistent with the recent results of cos- 
mological simulations which predict ste eper density profiles 
for more dynami c ally evolved clust ers (|Macci6 et al.ll20Qg| ; 
iDolag et al.ll2009l ; iDuffv et aLlboid ). 

The virial mass of the clusters correlates very well with 
the velocity dispersion, the X-ray temperature and the X-ray 
luminosity of the clusters. The slopes of the mass scaling re- 
lations agree with the estimates found in the literature and 
in the case of the velocity dispersion and X-ray tempera- 
ture are consistent with the prediction from the virial theo- 
rem. We used constraints on the virial mass and the concen- 
tration parameter to test the most up-to-date estimates of 
the erg parameter against the observational relation between 
the virial mass and the concentration parameter. Compar- 
ing the normalization of the mass-concentration relation for 
our cluster sample with the calibration from cosmological 
simulations we found tha t our data favours a rather high 
value, a s = 0.9±g;gg . Like iBuote et al.l (|2007l ). we found an 
excellent agreement with cr& from the 1st data release of the 
WMAP satellite. 

Our model of the projected phase-space density allowed 
us to constrain the profile of the orbital anisotropy in clus- 
ters. We showed at high significance level that galaxy orbits 
are isotropic in the cluster centres and radially anisotropic at 
the virial sphere, with /3 « 0.65 (a r /ag ~ 1.7). This finding 
is con sistent with the results obtained bv lBiviano & Katgert! 
(2004) for late type galaxies in nearby galaxy clusters and 
generalizes a number of estimates of th e global anisotropy 



of galaxy orbits in clusters (see e.g. Lokas et al 



isotropy 
1. 2006; 



IWoitak fc Lokas!l2007l ; lvan der Marel et al-lbOOoT ) 

Our constraints on the asymptotic values of the 
anisotropy parameter implies a growth of the anisotropy 
from f3 ~ at the cluster centre to /3 ~ 0.4 at the scale 
radius. Interestingly, this result agrees fairly well with re- 
cent esti mates of the ani s otropy for dark matter in galaxy 
clusters l|Host et al.ll2009l ; lHansen fc Piffarettill2007l l. 

The orbital anisotropy at the virial sphere of the low 
concentration clusters tends to be more radially biased than 
for the clusters of high concentration parameter. We suggest 
that this tendency is a consequence of a change in the galaxy 
orbits from radially biased at the early stage of cluster evo- 
lution to more isotropic for more dynamically evolved clus- 
ters. This scenario is supported by the studies of the orbital 
anisotropy for different galaxy populations which reveals a 
similar dichotomy between early and late type galaxies rep- 
resenting m ore and less relaxed gala xy populations in the 
clusters fsee lBiviano fc Katgertir2004h . 
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